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INTRODUCTION 


This  program  was  undertaken  as  a  portion  of  an  overall  effort  to  investigate  the  effects 
of  ballistic  damage  on  a  variety  of  metal  and  composite  rotor  blades.  When  a  rotor 
blade  spar  is  subjected  to  ballistic  impact,  the  resulting  damage  is  generally  an  irregularly 
shaped  hole  with  many  sharp  edges  that  produce  stress  concentrations.  Normally,  one 
sharp  edge  will  produce  the  most  critical  stress  concentration.  It  is  from  this  point  that 
the  first  fatigue  crack  will  initiate  and  propagate.  The  objective  of  this  program  was  to 
investigate,  through  analysis  and  test,  the  rate  of  crack  propagation  in  a  metal  blade  S|)ar 
from  an  induced  defect.  The  rotor  blade  evaluated  in  this  program  was  a  Bell  Helicop 
ter  Company  (BHC)  540  main  rotor  blade.  A  3/4  inch  diameter  hole  with  a  saw  cut 
that  was  stop  drilled  to  produce  a  controlled  crack  front  was  introduced  into  the  top 
surface  of  the  blade  spar.  The  blade  was  fatigue  tested  to  initiate  a  crack.  Fatigue  test 
ing  was  continued  at  maximum  level  flight  loads  and  the  crack  growth  was  monitored. 
The  experimental  data  was  then  compared  with  analytical  predictions  to  obtain  a  measure 
of  the  ability  to  predict  the  crack  growth  characteristics. 
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TFST  SPECIMEN  CONFIGURATION 


Th«  rotor  blade  testetl  was  a  Bell  Helicopter  Comjwny  (BHC)  Model  540  rnain  rotor 
blaile  artd  was  residual  protMuty  from  the  Morlel  540  Rotor  Blarie  Fatigue  Test  Program 
IRefererwre  1>.  The  blatie  is  approximately  19  feet  lortg  with  a  27  inch  chord  and  a 
9  1/3  percent  thick  symmetrical  airfoil  section.  Figure  1  shows  a  typical  cr^  swtion 
of  the  540  rotor  blade.  The  $»)ar  is  made  from  a  C  section  extrusion  o1  2024  aluminum 
alloy  amf  forms  the  basic  profile  of  the  leading  etige.  The  spar  $t>acer  is  alsii  a 
aluminum  alloy  extrusion  and  when  bonder!  in  place  controls  the  airfoil  thickness,  yre 
blatfe  from  the  s^rar  aft  is  composeif  of  aluminum  honeycomb  core  arnl  skirrs  Ironded  to 
the  siMir.  The  trailing  edge  is  a  strip  of  extruder!  aluminum.  The  sjrar  and  trailing  wlgr* 
strip  are  tatrererl  lietween  Statiom  BO  and  140  (see  Figure  2|  and  are  of  constant  chord 
outtHMrrf  of  Station  140. 

The  tip  end  of  the  blarfe  was  modified  by  BHC  for  testing  purt>oses.  A  22  iiKh  section 
of  the  tip  was  removerf  and  a  steiH»«d  aluminum  tfoubler  was  installeri  (s»h>  Figure  21. 

Tile  blade  test  stiecimen  had  been  ^>reviously  usihI  for  the  fatigue  test  evaluation  and 
investigation  of  the  H'«cer.  which  had  jiartially  rfelKiruferf  from  the  Sfiar  during  that 
testing'  Station  154  was  selected  as  the  test  section  for  this  test  tiecausr*  the  delwnrf  o1 
the  spacer  had  prot>agated  at  least  12  iiKhes  inlxiard  and  outlxsard  of  Station  154.  This 
precludtxf  protiagatiori  of  the  debond  through  the  test  section,  thus  eliminating  a  poten 
tial  unknown  variable  in  the  test  conditions. 
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Figure  1.  Cross  section  of  540  rotor  blade. 


'Gustafson  A.  J.,  and  Calapodas.  N.  J.,  Model  .^4(1  Rotor  HUidc  hatijtuc  7<Mf. 
USAAMRDL  Technical  Note  22.  Eustis  Directorate.  U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory.  Fort  Eustis.  Virginia.  January  1976,  AD  A021472. 
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Figure  2.  Modified  540  rotor  blade  for  testing. 


A  3/4-inch-diameter  hole  was  counterbored  in  the  top  side  of  the  blade  at  Station  154. 
A  chordwise  saw  cut  0.10  inch  in  length  was  made  from  the  hole  toward  the  leading 
edge  of  the  blade.  A  1 /64-inch  hole  was  drilled  at  the  tip  of  the  saw  cut  to  produce 
a  controlled  crack  tip  front.  The  geometry  of  the  induced  defect  in  the  blade  is  shown 
in  Figure  3. 


Figure  3.  540  rotor  blade  spar  geometry  at  Station  154. 
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TEST  EQUIPMENT 


The  Applied  Technology  Laboratory  rotor  blade  fatigue  test  machine  was  used  to  fabgue 
test  the  rotor  blade  specimen.  A  schematic  of  the  test  machine  is  shown  cn  Figure  4. 

The  centrifugal  force  is  applied  through  a  spring  system  and  the  flapwtse  bending  moment 
is  applied  through  a  load  actuator  as  shown.  The  rotor  blade  was  installed  in  the  test 
fixture  at  a  pitch  angle  to  minimize  the  chordwise  bending  moments. 


Control  of  the  fatigue  test  was  accomplished  through  a  closed  loop  servo  control  system 
computer  program  that  monitored  and  maintained  the  blade  bending  moment.  The 
resonant  frequency  was  found  by  manually  adjusting  the  frequency  up  and  down  in 
0.005  Hz  increments  and  was  then  maintained  automatically. 

The  blade  test  specimen  used  the  same  instrumentation  as  in  the  previous  testing,  as 
discussed  in  Reference  1.  The  moments  and  strains  developed  in  the  blade  were  wnsed 
by  the  strain  gages  and  were  recorded.  Data  acquisition  was  initiated  systematically  by 
the  control  system  computer  at  the  beginning  of  each  10,000  cycle  block  of  constant 
amplitude  load  cycles. 

After  completion  of  each  segment  of  constant  amplitude  load  cycles,  the  spar  was 
inspected  in  the  region  of  the  induced  defect  both  visually  and  with  ultrasonic  shear- 
wave  transmission  to  determine  the  location  of  the  crack  tip  front. 


8 


frame 


Figure  4.  Fatigue  test  machine-load  application. 


TEST  LOAD  DETERMINATION 


The  axial  and  bending  stresses  in  the  blade  at  Station  154  were  determined  from  the 
blade  stiffness  and  loads  data  contained  in  References  2  and  3,  respectively,  for  high¬ 
speed  level  flight.  The  blade  section  properties  and  high-speed  level  flight  loads  at 
Station  154  are  shown  in  Tables  1  and  2,  respectively.  A  stress  ratio,  R  =  0.74,  was 
calculated  for  this  flight  condition. 


At  Station  154  the  primary  steady  stress  is  due  to  the  centrifugal  force  and  the  pri¬ 
mary  alternating  stress  is  due  to  flapwise  bending.  Therefore,  to  simplify  the  test 
loading  condition  and  crack  propagation  predictions,  only  centrifugal  force  and  alterna¬ 
ting  flapwise  bending  loads  were  applied.  At  Station  154  the  centrifugal  force  is 
76,000  pounds.  To  maintain  the  stress  ratio,  R  =  0.74,  the  required  alternating  flap- 
wise  bending  moment  was  ±  7500  inch-pounds. 


TABLE  1.  BLADE  STIFFNESS  PROPERTIES 
AT  STATION  154 

_ STIFFNESS _ 

Flapwise,  35  x  10*  lb-in.  ^ 

Chordwsse.  3060  x  10*  lb-in.  ^ 

Torsion,  33.8  x  10*  lb-in. ^ 

Axial,  42.5  x  10*  lb 


TABLE  2.  HIGH-SPEED  LEVEL  FLIGHT  LOADS 
AT  STATION  154 

_ LOAD _ 

Flapwise  bending,  -3  +  7.7  x  10^  in.-lb 
Chordwise  bending,  88  +  26  x  10^  in.-lb 
Torsional  moment.  2.2  +  5.5  x  10^  in.-lb 
Centrifugal  force,  76  x  10^  lb 


^Structural  Analysis  of  540-01 1-] 00-13  Hub  and  Blade  Assembly  for  the  Model  209IAH- 
IJ,  UH-IC,  UH-IE,  AH-IG,  UH-IL,  UH-IM,  HH-IK  Helicopters,  Report  209-099-174 
Bell  Helicopter  Company,  Fort  Worth,  Texas,  November  1969. 

^Qualification  Load  Level  Survey  for  Improved  Main  Rotor  for  AH-IG,  Report  209-099- 
305,  Bell  Helicopter  Company,  Fort  Worth,  Texas,  June  1970. 
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CRACK  GROWTH  PREDICTIONS 


The  crack  growth  prectictions  were  made  using  Forman's  Equation  (Reference  4)  where 
the  crack  propagation  rate  is  a  function  of  the  stress  intensity  factor  range,  and  the 
stress  ratio,  R,  as  follows: 


da  ^  C  (J^K)” 

dN  '  (l-R)K^-JVK  (1) 

In  general,  the  stress  intensity  factor,  K,  is  a  function  of  the  maximum  stress  level,  o, 
and  the  half  crack  length,  a  ; 

K  =  (I  ^ '  ir  a 

Likewise,  the  stress  intensity  factor  range,  AK,  is 

AK  =  Au  ff  a 

wfiere  "  (l-RKr 


(2) 

(3) 

(4) 


For  2024  aluminum,  a  representative  value  for  n  is  3.0.  Substituting  Equation  (3)  and 
Equation  (4)  into  Equation  (1)  and  integrating,  the  following  equation  results: 


N, 


N,. 


Letting  Kf  =  Kc, 


Nc-No 


(5) 


16) 


The  above  equations  apply  to  the  case  of  the  centrally  cracked  infinite  plate  as  shown 
in  Figure  5.  However,  the  crack  growth  under  consideration  is  that  of  a  crack  propagat 
mg  from  one  side  of  a  hole  located  eccentrically  in  a  finite  width  plate.  To  treat  this 
case,  motfifications  to  the  theoretical  stress  intensity  factors  were  made  using  the 
methodology  of  Figgc  and  Newman  (Reference  5). 


■•Forman,  R.  0.,  Kearney,  V.  E.,  and  Engle,  R.  M.,  SunnrUiil  Amilysis  of  Crack  Propaga' 
(ion  in  (  vclic-l.iKiJcJ  Stnu  turcs,  Journal  of  Basic  Engineering,  September  1967. 

’Figge,  I.  E.,and  Newman,  J.  C.,  Jr.,  I'atifUic  Propayation  in  Srmcfnrcs  With 

SirniilatcJ  Kirct  I'orccs.  ASTM  Special  Technical  Publication  No.  415,  American  Society 
for  Testing  and  Materials,  Philadelphia,  Pennsylvania,  September  1967. 


In  this  case,  the  combination  of  crack  and  hole  is  assumed  to  tie  represented  by  an 
equivalent  crack  eccentrically  located  in  an  infinite  plate.  This  is  accomplished  by 
equating  the  stress  intensity  factors  for  an  infinite  plate  containing  a  hole  to  an  infinite 
plate  with  a  crack.  From  Reference  6,  these  stress  intensity  factors  are,  resfiectively: 

K  -  0  v/>r  /  fC^’/r)  (7) 

K  -  «  v''ir  ag  (8) 

where  the  function  fUVr)  shown  in  Figure  6  describes  the  influence  of  the  hole  on  the 
stress  intensity  factor.  Equating  these  and  solving  for  a^  results  in 

ag-nf(^/r)|'  (9) 

The  resulting  eiiuivalent  crack  length  is  found  to  be  nearly  equal  to  1/2  the  sum  of  the 
actual  crack  length  and  the  diameter  of  the  hole. 

The  stress  intensity  factor  can  be  further  modified  to  account  for  an  eccentrically  located 
crack  in  a  finite  width  plate  as  follows: 

K  »  0  ag  F  (e/b,  ag/b^)  (10) 

where  the  function,  F  (e/b,  ag/b^),  shown  in  Figure  7,  is  the  boundary  correction 
factor  for  the  eccentrically  locat^  crack  in  a  finite  width  plate.  If  Equation  (10)  is 
then  equated  with 

K  “  ()  y/jr  Bg  (11) 

where  Sg  is  the  effective  half  crack  length  (taking  into  account  the  presence  of  the  hole, 
eccentricity,  and  finite  width)  then 

Sp  •  nMt’/r)  FMe/b,  ag/b,)  (12) 

The  rotor  blade  test  specimen  geometry  was  idealized  as  shown  in  Figure  8.  Appendix 

A  shows  the  calculated  effective  half-crack  length  and  resulting  crack  growth  rate  and 
cycles  to  failure  for  initial  crack  lengths,  f  o,  varying  from  0.25  and  2.75  inches.  This 
calculation  uses  Equations  (1)  and  (6)  restHK;tivelv,  together  with  the  stress  intensity 
factor  of  Equation  (1).  The  stress  conditions  used,  which  correspond  to  a  centrifugal 
force  of  76,000  pounds  and  an  alternating  flapwise  bending  moment  of  ±  7500  inch 
pounds,  are 


o  -  21403  Ib/in.^ 
A<i  -  5566  Ib/in.^ 


•Paris,  P.  C.,and  Sih,  G.C.M.,  Sin'ss  Atutlvsis  of  Cnu  ks,  ASTM  Sgrcial  Technical 
Publication  No.  381,  American  Society  for  Testing  and  Materials,  Philadelphia,  Pennsyl 
vanla,  June  1d64. 
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Th«  values  of  the  constants  and  C  determined  by  Hudson  (Reference  7)  for  2024 
aluminum  were  selected  as  representative  values  for  the  blade  specimen.  They  are 

Kg  -  100320  Ib/in.''- 
C  -  3.22  X  10' 


Figure  6.  Stress-intensity-factor  coefficient  for  a  single  crack 
emanating  from  a  circular  hole  in  a  uniaxial  stress 
state. 

^Hudson,  C.  M.,  Effect  of  Stress  Ratio  on  Fatigue-Crack  (Irowth  in  7()75-T6  ami 
J024‘Tt>  Aluminum-Alloy  Specimens,  NASA  Technical  Note  0-5390,  National  Aeronautics 
and  Space  Administration,  Washington,  D.  C.,  August  1969. 
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Figure  7.  Correction  factor  for  eccentricity  and  finite  width. 


TEST  PROCEDURE 


The  test  specimen  was  installed  in  the  rotor  blade  test  fixture  and  a  steady  centrifugal 
force  of  76,000  pounds  and  an  alternating  flapwise  moment  of  ±  7500  inch-pounds  at 
Station  154  was  applied.  The  specimen  was  loaded  in  blocks  of  10,000  cycles  each. 

Since  a  crack  had  not  initiated  after  three  load  blocks,  the  alternating  flapwise  bending 
moment  was  increased  to  ±  11,000  inch-pounds  at  Station  154  in  an  attempt  to  initiate 
a  failure.  During  this  load  block  a  crack  initiated  and  propagated  a  total  of  0.15  inch. 
The  alternating  flapwise  moment  was  then  reduced  to  ±  7500  inch-pounds  and  main¬ 
tained  at  that  level  during  each  subsequent  block  of  10,000  load  cycles. 

The  specimen  was  inspected  in  the  area  of  the  induced  effect  visually  on  the  outside 
and  on  the  inside  through  a  mirror  system.  The  crack  tip  front  was  located  using  shear 
wave  ultrasonics  and  dye  penetrant.  Once  the  crack  tip  front  was  located,  the  length  of 
the  crack  was  then  measured. 

The  inspection  of  the  spar  was  conducted  after  each  block  of  load  cycles  through  block 
16  and  then  after  every  other  block  of  load  cycles  through  block  36.  The  final  inspec¬ 
tion  prior  to  ultimate  failure  was  made  after  block  42. 


TEST  RESULTS 


The  load  history  and  measured  crack  length  for  the  test  specimen  are  set  forth  in  Appen 
dix  B.  The  relationship  between  the  actual  crack  length,  t  ,  and  the  number  of  applied 
load  cycles  is  shown  in  Figure  9,  where  the  initial  condition,  =  O,  is  assumed  to  be 
the  last  cycle  of  load  block  4.  The  corresponding  initial  crack  length,  (  o,  is  then  0.25 
inch. 

Using  the  methodology  of  the  previous  selections,  the  effective  half-crack  length,  a^, 
was  determined  for  each  crack  length,  t.  The  resulting  relationship  between  ag  and  the 
number  of  applied  load  cycles  is  shown  in  Figure  10. 

The  crack  propagation  rate  in  the  spar  was  determined  using  the  following  methods: 

1.  A  best-fit  curve  was  drawn  through  the  data  shown  in  Figure  10.  The  slope 
of  the  curve  was  then  determined  using  finite  difference  techniques.  Figure 

11  shows  the  resulting  relationship  between  the  crack  growth  rate,  dag/dN,  and 
the  stress-intensity-factor  range,  AK.  For  comparison  the  crack  growth  rate, 
dig/dN,  determined  using  Forman's  Equation  (Equation  10)  is  also  shown  in 
Figure  11. 

2.  A  best  fit  curve  was  drawn  through  the  data  shown  in  Figure  9.  The  slope  of 
the  curve  was  determined  using  finite  difference  techniques.  Figure  12  shows 
the  resulting  relationship  between  the  crack  growth  rate,  d  f  /disi,  and  the 
stress  intensity-factor  range,  AK. 

3.  The  failed  spar  was  sectioned  as  shown  in  Figure  13.  The  fracture  surface  was 
examined  at  the  locations  shown  using  the  electron  scanning  microscope  at 
NASA/Langley  Research  Center.  Appendix  C  shows  a  series  of  fractographs  of 
the  areas  of  the  fracture  surface  that  was  examined. 

From  the  fractographs  the  fatigue  striations  were  counted  and  the  crack  growth  rate 
d£/dN  was  subsequently  determined.  The  resulting  crack  growth  rates  are  shown  in 
Figure  12. 

In  each  case  above,  the  stress-intensity-factor  range  is  based  on  the  effective  half-crack 
length,  ig.  Examination  of  the  crack  growth  rate  curves  shown  in  Figure  1 1  shows  a 
reasonable  correlation  between  the  crack  growth  rate  interpolated  from  the  test  data  and 
the  rate  predicted  from  Forman's  Equation. 

Examination  of  the  data  shown  in  Figure  12  reveals  a  similar  trend  between  the  crack 
growth  rates  determined  from  the  slope  of  the  best-fit  curve  of  Figure  9  and  the  crack 
growth  rate  determined  from  the  striation  count.  The  distinctive  characteristic  of  the 
data  shown  in  Figure  12  is  that  the  crack  growth  rate  is  initially  high  and  then  drops 
off  before  continuing  to  increase.  Recalling  that  the  crack  growth  was  originally  initiated 
at  the  higher  load  level  of  ±  11,(X)0  inch-pounds  and  the  load  was  then  reduced  to 
±  7500  inch-pounds,  it  appears  that  the  phenomenon  referred  to  as  "delay  retardation"  in 


Reference  8  occurred.  That  is  to  say  that  the  crack  growth  rate  decreased  after  the 
higher  load/overload  cycling.  At  a  point  where  the  crack  length,  had  reached  approxi¬ 
mately  1  inch,  the  crack  growth  had  resumed  a  normal  crack  growth  trend. 

The  "delay  retardation"  effect  does  not  appear  in  the  crack  growth  rate  versus  stress- 
intensity-factor  range  shown  in  Figure  11.  This  may  be  a  result  of  the  curve  smoothing 
procedure.  The  lower  initial  crack  growth  rate  compared  to  that  predicted  by  Forman's 
Equation,  however,  may  be  due  to  the  "delay  retardation"  effects.  The  difference  in 
the  slope  of  the  curves  shown  in  Figure  10  may  be  attributed  to  variations  in  the 
selected  material  constants  C  and  n. 

The  critical  crack  length,  I  in  testing  was  found  to  be  4.75  inches  based  on  visual 
examination  of  the  fracture  surface  after  the  test  was  completed.  In  comparison,  the 
predicted  critical  crack  length  using  Forman's  Equation  was  4.85  inches.  The  fact  that 
the  actual  critical  length  was  slightly  less  than  predicted  would  seem  to  indicate  that  the 
actual  critical  stress-intensity  factor  is  somewhat  less  than  the  representative  value  selected. 


*  Hertzberg,  R.  W.,  Deformation  ami  Fracture  Mechanics  of  Fugineerinf:  Materials,  New 
York,  John  Wiley  and  Sons,  1976. 


Figure  12.  Crack  propagation  rate  versus  stress  intensity-factor  range. 
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2. 


CONCLUSIONS 


Based  on  the  data  obtained  from  the  testing  of  this  single  test  specimen,  it  is  concluded 
that: 


Forman's  Equation  can  be  used  to  predict  the  crack  growth  behavior  of  a 
controlled  defect  in  a  metal  blade  spar  if  the  appropriate  modifications  are 
made  to  the  stress-intensity  factor  to  incorporate  corrections  for  hole 
influence,  eccentricity,  and  finite  width. 
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APPENDIX  B 

TEST  LOAD  AND  CRACK  GROWTH  RECORD 


Nominal  Altarnata 


Load  Block 
Numbar 

Flap  Banding 
(in.-lb) 

Crack  Propagation  (in.) 
During  Block  Total 

1 

*  7,500 

0 

0 

2 

+  7300 

3 

+  7300 

4 

t  11.000 

0.15 

0.25 

S 

±  7300 

0.10 

0.35 

6 

♦  7,500 

0 

0.35 

7 

±  7,500 

0 

0.35 

8 

+  7300 

0.19 

0.54 

9 

±  7300 

0 

0.54 

10 

i  7,500 

0.03 

0.57 

11 

i  7300 

0.03 

0.60 

12 

i  7,500 

0 

0.60 

13 

i  7300 

0.18 

0.76 

14 

±  7,500 

0.76 

15 

±  7,500 

0.06 

032 

16 

i  7,500 

0.06 

038 

17,18 

+  7,500 

0.03 

0.91 

19,20 

±  7,500 

0.10 

1.01 

21,22 

+  7300 

0.06 

1.07 

23,24 

±  7300 

0.12 

1.19 

25,26 

±  7300 

0.10 

1.29 

27,28 

i  7300 

0.10 

1.40 

29,30 

t  7300 

0.12 

1.52 

3M2 

t  7300 

0.12 

1.64 

Commenti 


Tatt  load  increawd  to  initiata  crack. 
Crack  initiated,  running  to  leading  edge. 
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Load  Block 
Number 


Nominal  Altamata 

Flip  Banding  Crack  Propagation  (in.) 

_ 6uring  Block  Total  _  Comment* 


33,34 

♦  7,500 

0.35 

1.89 

35,36 

*  7,500 

0.25 

2.14 

37 

39,40 

41,42 

i  7,500 

0.71 

2.85 

43 

±  7,500 

1.90 

4.75 

Continued  running  •  no  impaction. 


Ultimata  failure. 


NOTES; 

1.  Steady  axial  load  was  nominally  76,000  pound*. 

2.  The  load  block*  con*i*t*d  of  10,000  cycle*  each. 

3.  The  crack  growth  wa*  determined  u*ing  ultraionic  (hear  wave  impaction  technique*. 

crack  growth  wa*  determined  after  each  block  or  group  of  block*  a*  ihown  above. 

5.  The  alternating  bending  wa*  removed  and  the  axial  load  reduced  to  30,000  pound*  for  each  impaction. 

6.  Total  crack  propagation  length  ihown  include*  0.10-inch  initial  *aw  cut. 
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APPENDIX  C 

SCANNING  ELECTRON  FRACTOGRAPHS 
OF  FRACTURE  SURFACE  OF  540  ROTOR  BLADE  SPAR 


NOTE:  Th«  l«t  mo  thr»*<«9it  numb*rt  on  Mch  frictogrtph  ar*  the  Mction  and  location  number*  identified 
in  Figure  10. 
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LIST  OF  SYMBOLS 


a 

= 

half-crack  length 

= 

effective  half-crack  length 

h 

= 

effective  half-crack  length 

f(  > 

= 

a  function  of 

= 

actual  crack  length 

1 

c 

= 

critical  actual  crack  length 

o 

= 

initial  actual  crack  length 

n 

= 

numerical  exponent 

r 

= 

radius 

C 

= 

material  constant 

K 

= 

fracture  mechanics  stress-intensity  factor 

AK 

= 

stress-intensity-factor  range 

Xc 

= 

critical  stress  intensity  factor  for  fracture 

N 

cycle  number 

No 

= 

initial  cycle  number 

Nc 

= 

critical  cycle  number  at  crack  instability 

R 

= 

ratio,  minimum  stress  to  maximum  stress 

da/dN 

= 

crack  extension  per  cycle  of  load 

a 

= 

applied  stress  normal  to  crack 

Aa 

applied  stress  range 
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